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Euclidean Maxwell theory in the presence of boundaries. II 
Abstract. C-function regularization is applied to complete a recent analysis of the quan- 
tized electromagnetic field in the presence of boundaries. The quantum theory is studied 
by setting to zero on the boundary the magnetic field, the gauge-averaging functional and 
hence the Faddeev-Popov ghost field. Electric boundary conditions are also studied. On 
considering two gauge functionals which involve covariant derivatives of the 4-vector po- 
tential, a series of detailed calculations shows that, in the case of flat Euclidean 4-space 
bounded by two concentric 3-spheres, one-loop quantum amplitudes are gauge independent 
and their mode-by-mode evaluation agrees with the covariant formulae for such amplitudes 
and coincides for magnetic or electric boundary conditions. By contrast, if a single 3-sphere 
boundary is studied, one finds some inconsistencies, i.e. gauge dependence of the ampli- 
tudes. 
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Euclidean Maxwell theory in the presence of boundaries. II 
1. Introduction 

Despite the lack of a mathematically consistent theory of quantum gravity, the elliptic 
boundary-value problems occurring in quantum cosmology have recently shed new light 
on the whole quantization programme for gauge fields and gravitation in the presence 
of boundaries [1-4]. Boundary effects play a key role in the path-integral approach to 
quantum gravity [1] , in the problem of boundary conditions for the quantum state of the 
universe [5], and in comparing different quantization and regularization techniques for 
gauge fields and gravitation [1-2]. The latter problem provides the main motivation for 
our paper. In fact many efforts have been produced in the literature to understand the 
relation between canonical and manifestly gauge-invariant approaches to quantum field 
theories, as well as to compare mode-by-mode and covariant formulae for the evaluation 
of quantum amplitudes. 

In particular, in some papers dealing with the calculation of the scaling factor of the 
wave function of the universe in quantum cosmology, discrepancies were found between 
results obtained by covariant and non-covariant methods. It is well-known that this scaling 
factor coincides with the Schwinger-De Witt coefficient A2 in the heat-kernel expansion 
[6]. Moreover, this factor can be calculated by using the generalized Riemann ^-function 
technique [1]. Within this framework, the prefactor is expressed through the C(0) value, 
while C'(0) yields the full expression for the one-loop effective action. It was noticed 
that, for fields with non-zero spin, calculations of the Schwinger-De Witt coefficient A2 [6] 
by using general covariant formulae for Riemannian 4-manifolds with boundaries [7] give 
results which differ from those obtained by using ^-function technique when one restricts 
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the theory to its physical degrees of freedom [1]. An analogous phenomenon was noticed 
in [8], where C(0) was calculated for gravitons on the full Riemannian de Sitter sphere. In 
[9] the hypothesis was put forward that the reason of discrepancies lies in the impossibility 
to perform a 3+1 decomposition on the Riemannian 4-manifolds under consideration. 
Indeed, this decomposition is necessary for the separation of physical degrees of freedom. 
Thus, the direct calculation of C(0) in terms of physical degrees of freedom seems to be 
inconsistent. The ^0) calculation for fermionic fields was then carried out in [9] on the 
part of flat Euclidean 4-space bounded by two concentric 3-spheres. It was shown that the 
discrepancy disappears in this case. 

However, to understand discrepancies for the electromagnetic field and other gauge 
theories, it is necessary to take into account ghost modes and non-physical degrees of 
freedom, whose contributions may survive in a non-trivial background even in the unitary 
gauges. We here restrict ourselves to a mode-by-mode analysis of the scaling factor in 
relativistic gauges about fiat Euclidean 4-space with one or two 3-sphere boundaries. One 
then faces the following problems, here described in the case of Euclidean Maxwell theory, 
which is the object of our investigation. 

(i) Choice of background 4-geometry. This can be fiat Euclidean 4-space, or a curved 
Riemannian manifold providing the index of the Dirac operator vanishes and no further 
obstructions to having a unique, smooth solution of the classical, elliptic boundary-value 
problem can be found. [Knowledge of the index of the Dirac operator ensures one under- 
stands what happens for second-order elliptic operators as well] 
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(ii) Choice of boundary 3-geometry. Motivated by quantum cosmology, this is taken 
to be a 3-sphere, or two concentric 3-spheres. These choices are necessary to have a unique 
smooth solution of the corresponding classical boundary-value problem for fields of various 
spins, and to avoid singularities at the origin of the background 4- geometry (see below). 

(iii) Choice of boundary conditions. They can be magnetic, which implies setting 
to zero on the boundary the magnetic field, the gauge-averaging functional and hence 
the Faddeev-Popov ghost field. They can also be electric, hence setting to zero on the 
boundary the electric field, and leading to Neumann conditions on the ghost [1] . 

(iv) Choice of gauge-averaging functional. Here we focus on the gauge- averaging func- 
tional first proposed in [1-2]: ^e{A) = ^^"^V^Ah—Aq Tr K {K being the extrinsic-curvature 
tensor of the boundary), and on the Lorentz functional $l(^) = ^^^V'^^^, as a first check 
of gauge independence of quantum amplitudes in a mode-by-mode analysis of the quantized 
electromagnetic field. 

The mode-by-mode analysis is performed by relying on the familiar expansions of the 
components of the 4- vector potential on a family of 3-spheres centred on the origin [1-3] , 
i.e. 



Ao{x,T)^Y.^r^{r)Q^''\x) 



(1.1) 



n=l 



oo 

Ak{x,T) = J2 [fn{r)Sl!'\x)+gn{T)Pt''\x)] for aU /c = 1,2,3 (1.2) 



n=2 



where Q^"'\x) , sj^\x) , P^"'\x) are scalar, transverse and longitudinal vector harmonics 
on respectively. Note that, however, normal and tangential components of A/j, are 
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only well-defined at the 3-spliere boundary where t = t^, since a unit normal vector 
field inside matching the normal to at the boundary is ill-defined at the origin of the 
coordinate system for fiat Euclidean 4-space [9]. Hence the geometrical meaning of (1.1)- 
(1.2) as normal and tangential components of the 4- vector potential inside the 3-sphere 
boundary remains unclear, unless one studies an elliptic boundary-value problem where 
fiat Euclidean 4-space is bounded by two concentric 3-spheres of radii t+ and t_ , say. The 
results of our calculations show that a proper study of non-physical degrees of freedom and 
ghosts (which do not compensate each other), together with the consideration of manifolds 
possessing a consistent 3-|-l decomposition, eliminates the discrepancies between covariant 
and non-covariant formalisms. 

Our paper is organized as follows. Section 2 studies one-loop amplitudes by choosing 
the gauge-averaging functional $£(^4) (see above), and imposing magnetic (or electric) 
boundary conditions on 3-spheres. The one-boundary and two-boundary problems are 
analyzed. Section 3 repeats the investigation of section 2 in the case of the gauge functional 
^l{A). Results and open problems are presented in section 4. Details relevant for ^(0) 
calculations are described in the appendix. 

2. One-boundary and two-boundary problems in the Esposito gauge 

In this section we first evaluate ({0) for the electromagnetic field on the fiat 4-dimensional 
Euclidean background bounded by a 3-sphere. We choose the magnetic boundary condi- 
tions described in the introduction and carry out our calculations in the Esposito gauge 
[1-2] . For this purpose, we begin by studying the coupled eigenvalue equations for normal 
and longitudinal components of the electromagnetic field obtained in [1-2] . They have the 
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form (hereafter we set to 1 the parameter appearing in the Faddeev-Popov action (2.3) of 

[2]) 

An 9n{r) + Bn Rn{r) ^ (2.1a) 
dngnir) + DnRn{r) = Q (2.16) 

where 



^ Id (n2-l) ^ ^^^^ 

dr^ T dr 



ij„ . -i!^ (2.3) 



C„ = — L (2.4) 



5 - + ^ ^ ^^'-^) I A f25) 



To study non-physical modes it is convenient to diagonalize the operator matrix 



An Bji 
Cn Dn 



Hence we fook for a diagonalized matrix in the form 

4-^(4.) "f')- 

The matrix '^'i^ creates the hnear combinations of functions RniT) and Qni^) 

which can be found from decoupled equations, whilst the matrix ( ^ selects these 



decoupled equations. 
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Setting to zero the off-diagonal matrix elements of (2.6), and defining v = +-y/n^ — |, 
one finds equations solved by Vn = —an, Wn = —Pn, where 



an{T)={-l±-'] T (2.7) 



Choosing the pair of solutions with upper sign for a^, and lower sign for /3n [the opposite 
choice of signs gives the equivalent system of operators whilst the choice of coinciding signs 
for an and (3n implies the degenerate system of equations] one finds general basis functions 
for Rn{T) and Qni^) in the form 

Qnir) = Ci /._i/2(VAt) + C2(v- 1/2) /.+i/2(VAt) (2.9) 

Rn{T) = ^ (C, (^7^^ /.-1/2(VAt) + C2 /.+1/2(VAt)^ (2.10) 

where Ci and C2 are constants. 

In our particular gauge, magnetic boundary conditions imply Dirichlet boundary con- 
ditions for gni^) and Neumann boundary conditions for Rni^) [1-2]. The resulting equa- 
tions lead to a 2 X 2 matrix, hereafter denoted by X, whose determinant has to vanish 
to find non-trivial solutions. Such an eigenvalue condition for normal and longitudinal 
components of the electromagnetic potential is best studied by using the algorithm of [4] . 
It is known that C(0) can be expressed as 

C(0) = /log + /pole(oo) - /pole(O) (2.11) 



Euclidean Maxwell theory in the presence of boundaries. II 



where [4] 



oo 



/pole(M2) 



I{M^,s)= d{n) 



n-^' log UM^) 



s 



+ /«(M2) + 0(s). 



(2.12) 



With our notation, d{n) is the degeneracy of the eigenvalues parametrized by the integer n, 
and fn{M'^) is the corresponding eigenvalue condition. Moreover, on analytic continuation, 
/log = I^g is the coeflBcient of log M from I{M'^, s) as M — > oo, and Ipo\e{M'^) is the residue 
at s = 0. The condition det 1 = should be studied after eliminating fake roots M = 0. 
To obtain that it is enough to divide det X by the minimal power of M occurring in the 
determinant. It is easy to see by using the series expansion for modified Bessel functions 
[10] that such a power is M'^^~^. 

We begin with the calculation of Iiog for normal and longitudinal modes of the elec- 
tromagnetic field together with ghosts. Using uniform asymptotic expansions for modified 
Bessel functions [10] we can see that the only terms in the logarithm of det 1 divided by 
M^^~^ which are proportional to logM have the form —21^ logM whilst the ghost eigen- 
value condition, divided by M^, gives analogous terms — (2z^ -|- 1) logM which contribute 
to /log with the opposite sign. Hence we can write 



where is the dimension of the irreducible representation for scalar hyperspherical har- 
monics. Summation in (2.13) is carried out by the method of (^-function regularization 
(see for details [4]). The infinite sum starts from n = 2 since effects of zero- modes for 




(2.13) 
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ghosts and normal photons should be calculated separately whilst the longitudinal (n = 1) 
photon is absent. In [1-2] it was found that the contribution of the decoupled normal mode 
with magnetic boundary conditions is 



C(0)decoupled mode — . • 

(2.14) 



It is easy to calculate the contribution to C(0) resulting from ghost zero-modes by substi- 
tuting into the corresponding expression for /log the value n = 1. One finds 

C(0)ghost zero— modes — 1- (2-15) 

Now we have to calculate /poie(oo) and /poie(O). As shown in the appendix, the 
structure of the term generating /poie(oo) is j^^iYJi)- Taking the logarithm of this expression 
and expanding it in inverse powers of n we can pick out the coeflBcient of the term 1/n in 

2 o 

the expression ^^ log as n — > oo, and find 



/pole(oo) = -^. (2.16) 



Now we can calculate Ipoie 

(0) simply by using the usual series expansion of Bessel functions 
[10] in the limit M — > 0. The calculation relies on the Stirling formula for the F-function 
appearing in the series expansion for Bessel functions [10] . Omitting the details we write the 
result for /poie(O) for normal and longitudinal modes of the electromagnetic field together 
with ghosts (see appendix) 

5 

-^pole (nonphys and ghosts) (0) (^■''"7) 
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Now combining the results (2.13)-(2.17), one has 



53 

C(0) normal, longitudinal, ghosts — gg- (2.18) 



Remarkably, the total contribution of non-physical degrees of freedom and ghosts does not 
vanish. Adding to (2.18) the contribution —j^ of physical modes obtained in [3] one finds 



We now evaluate C(0) for the electromagnetic field on the flat 4-dimensional Euclidean 
background bounded by two concentric 3-spheres. To begin with, let us calculate the 
contribution to C(0) of physical degrees of freedom. The basis functions for them are now 
the linear combination /n(T) = Ci IniMr) + C2 Kn{MT) which should vanish at the 
3-sphere boundaries of radii tj^ and t_ respectively, where T-f > t_. This leads to the 
eigenvalue condition 

/- - /+ i^- = (2.20) 

where /- = /„(Mt_),/+ = In{MT+),K- = Kn{MT-),K+ = Kn{MT+). Using series 
expansions for modified Bessel functions one can see that the eigenvalue condition (2.20) 
has no fake roots. Since in (2.20) the coeflicients of products of Bessel functions are 
independent of n, one finds that /poie(oo) = 0. Bearing in mind that the dimension of 
irreducible representations for transverse vector hyperspherical harmonics is 2 (n^ — 1) 
[3], where n = 2, . . ., and using the hypergeometric expansions for and K^, one can 
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easily show that /poie(O) = 0. Using the asymptotic expansions for In and one can also 
calculate /log as 

/iog = -E(^'-l) = -9 (2-21) 



n=2 



and correspondingly 



C(0) transversal photons — (2.22) 



It is not difficult to consider also the contribution of ghosts. In this case one has an 
eigenvalue condition coinciding with (2.20) if we use v = -\-\/v?- — 3/4 instead of n. Hence 
all 3 contributions to (^(0) vanish, which implies 

C(0)ghosts = 0. (2.23) 

The next problem is the calculation of the contribution to C(0) of the decoupled 
normal mode, which has the form [1-2] Rx{t) = Cx \ I\{Mt) + C2 ^ Ki{Mt). The 
derivative of this function should vanish at the 3-sphere boundaries. The determinant of 
the corresponding 2x2 matrix should vanish and this equality gives, as in the previous 
cases, the eigenvalue condition. Such a determinant has no fake roots. Thus, by using the 
uniform asymptotic expansions of Bessel functions one can see that the I\og value is — |. 
However, since our decoupled mode is non- vanishing for r e [t_, t+], one deals with a zero 
eigenvalue corresponding to a non-zero eigenfunction satisfying boundary conditions. The 
number Nd = 1 of such decoupled normal modes contributes to the full ('(O) value. Hence 
one finds 

C(0)decoupled mode = -^log + Nd ~ ~2 ^ ~ 2' (2.24) 
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We now evaluate the contribution of the coupled normal and longitudinal modes of 
the electromagnetic field. Since in the two-boundary case the singularity at the origin is 
avoided, both /- and K- functions contribute to gauge modes. Hence the general solutions 
for gn{T) and Rnij) are 



After substitution of (2.25)-(2.26) into magnetic boundary conditions at the 3-sphere 
boundaries one has a system of 4 equations. The determinant of the corresponding 4x4 
matrix should vanish. In such a determinant, the smallest power of M is M~^. Thus, to 
avoid the appearance of fake roots in the eigenvalue condition, it is necessary to multiply 
it by M^. Taking into account the formulae for the asymptotic expansions of Bessel func- 
tions one can see that the coefficient of log M in the logarithm of our eigenvalue condition 
vanishes, hence /log vanishes as well. Let us now calculate /poie(oo)- Just as in the previ- 
ous problem, this value is determined by n-dependent coefficients in the determinant and 

can be calculated from the expression ^ log ^ (1^+1/2)2 ) along the lines described in the 

appendix. One has /poie(oo) = - h. ^ 



9n{T) - Ci /,_i/2(Mt) +C2{v- 1/2) /,+i/2(Mt) 



+ CS i^.-l/2(MT) + C4 (iV - 1/2) K,+i/2(Mt) 



(2.25) 





(2.26) 



13 



Euclidean Maxwell theory in the presence of boundaries. II 
Now using the usual hypergeometric expansions for modified Bessel functions one 
obtains /poie(O) = and since the values of /poie at oo and at compensate each other, 
and /log = 0, one finds 

C(0)normal and longitudinal — 0- (2.27) 

Combining the results (2.22)-(2.24) and (2.27) one has 

C(0) = (2.28) 

which coincides with the covariant result. Of course, the covariant ^(0) value is zero, since 
the volume contribution on the flat background vanishes whilst two-boundary contributions 
compensate each other (look at the corresponding formulae in [7]). 

On imposing electric boundary conditions, which are motivated by supersymmetric 
quantum field theory [1] , the situation for coupled normal and longitudinal modes is just 
opposite to that in the magnetic case, since normal modes and the normal derivatives 
of longitudinal modes vanish at the 3-sphere boundaries [1-2]. Defining = (^j^^^, 
bi, = (v — 1/2), the corresponding eigenvalue condition is the vanishing of the determinant 
of the matrix 



/ 


a„If^ 


-^6^+1 






\ 






/+ 






(Ik 




&.(/,; + /,;+2) 










-l+^^+l) 








J 



This yields 



/log = /pole (oo) = /pole (0) = (2.29) 
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which imphes 



normal and longitudinal — 0- 



(2.30) 



Hence one finds 



C(o) = 0. 



(2.31) 



Such a C(0) value coincides with the covariant result (as in the case of magnetic boundary- 
conditions, the contributions of two 3-sphere boundaries in covariant formalism cancel each 
other) . 

It is easy to carry out the corresponding calculations in the Lorentz gauge (section 
3). They yield again (^(0) = 0. Thus, we have shown that in the case of electric boundary 
conditions (just as in the case of magnetic ones) on the flat Riemannian 4-manifold with 
two 3-sphere boundaries leading to a well-defined 3+1 split of the 4- vector potential, results 
of covariant and mode-by-mode formalisms coincide. 

3. One-boundary and two-boundary problems in the Lorentz gauge 

In the Lorentz gauge, coupled eigenvalue equations for normal and longitudinal modes are 
again described by a system of the kind (2.1a)-(2.1b). However, the four operators are 
replaced by other operators En, Fn,Gn, Hn respectively, where the operator En coincides 
with An written in section 2, whilst Gn and Hn are different and take the form 



2 (n^ - 1) 



(3.1a) 



T 




(3.16) 
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Hn=-r^ + -—-^ ^ + A„. (3.1c) 



dr^ T dr 



T 



The corresponding system of equations can be diagonalized in the same way as the anal- 
ogous system (2. la)- (2. lb). Hence one finds a couple of Bessel-type equations, giving the 
general solutions for gn{T) and Rn{j) in the one-boundary problem as 

gn{r) = Ci /n+i(MT) + C2 (n + 1) In-iiMr) (3.2) 

Rn{r) = ^(-Cl J^^^ In+l{MT) + C2 In-l{MT)). (3.3) 

In the magnetic case we are studying, one sets to zero on the boundary the gauge-averaging 
functional. If this is the Lorentz functional, one is choosing Dirichlet conditions for 
modes and Robin conditions for modes, i.e. gn{T+) = 0, ^^(t+) -|- :^i?n(T+) = [1-2]. 
Such a system leads to the eigenvalue condition 

/ It+i {n + l)I+_, \ 

y {n-l) \{Mt+/2) ^ ^ j \{Mr+/2)^^n-2^^n j j 

To get rid of fake roots it is necessary to divide the determinant (3.4) by M^"^"^, and then 
we can calculate /log as 

/,., = -f;„3 = H5. (3.6) 

n=2 

/poie(oo) can be obtained by extracting the n-dependent coefficients in the determinant 
(3.4), which gives \ log ^ {n-i) ) • ^Vom this expression one obtains (see appendix) 



-fpoie(oo) = ^. (3.6) 
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Now we can calculate /poie(O) by taking the logarithm of the determinant (3.4) in the 
limit M — > and expanding it in inverse powers of n. The result is (see appendix) 



. ^ 1 1 121 , „^ 

/pole = + - = . 3.7) 



Combining together the results (3.5)-(3.7) one finds 



37 

C(0)normal and longitudinal — (^•^) 



The equation for the decoupled normal mode implies Ri — ^ I2{Mt), and the mag- 
netic boundary condition leads to 



3 

C(0)decoupled = — ^- (3.9) 



The eigenvalue condition for ghosts in the Lorentz gauge coincides with the one for 
scalar fields and the corresponding (^(0) value can be obtained from the well-known result 
for a scalar field subject to Dirichlet boundary conditions [1]. Taking into account the 
change of sign due to the fermionic nature of ghosts corresponding to spin-1 fields [1-2] 
one has 

C(0)ghosts = ^. (3.10) 
Combining the results (3.8)- (3. 10) one obtains 



C(0)normal, longitudinal, ghosts — ' (^-H) 
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Correspondingly, the total C(0) value including the contribution —j^ of physical modes is 



(3.12) 



Remarkably, (3.12) agrees with the corrected C(0) value obtained in [7], which relies on the 
analysis by Vassilevich [11]. However, (3.12) differs from the mode- by- mode result (2.19). 
This discrepancy seem to originate from the ill-definiteness of the 3-1-1 decomposition of 
the 4- vector potential on the manifold under consideration [9] . 

In the two-boundary case, the eigenvalue condition for the coupled normal and longi- 
tudinal modes in the Lorentz gauge on the flat Riemannian 4-manifold with two concentric 
3-sphere boundaries is (on imposing magnetic boundary conditions) 



det 



V 



T+ 

1 ( "^+1 

(n-1) \ Mt_ /2 
1 

(n-1) V M-r_|_/2 



n+2 



) 



(n + l)/,_i 
( 



M-r_ /2 



Z+In) 



't-1 



+/+_2+-f^) 



i-i) ( 



IK 



Mt_ /2 



n n+2 / 

/ 2X+ , , 
(n-1) V M-r_|_ /2 

n+2 / 



(n + l)K-_i\ 
(n+l)K+_i 



IK 



\ Mt_ /2 

-KZ ^-K 



-A 



Mt_|_ /2 



= 0. 



-Kt) I 



(3.13) 



The contribution to /log of the determinant (3.13) is 0. The function determining the 
behaviour of /poie(-^^) as M — >^ oo and n — > cxo has the form ^ log ^^^^^jp-j and corre- 
spondingly /poie(oo) = \. After taking the logarithm of the determinant (3.13) at M = 
and expanding it in inverse powers of n one obtains /poie(O) = \- Thus, one finds again 
compensation of /poie(oo) and /poie(O) and C(0)normai and longitudinal = 0. Finally one obtains 
the total result 



C(o) = 



(3.14) 
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coinciding with that obtained in section 2 and with the covariant one. 

4. Concluding remarks 

The main results of our investigation are as follows. 

First, we have proved that, in the case of flat Euclidean 4-space bounded by two 
concentric 3-spheres, the mode-by-mode analysis of one-loop quantum amplitudes for Eu- 
clidean Maxwell theory agrees with the covariant formulae used in [7]. Moreover, such 
Faddeev-Popov quantum amplitudes are indeed gauge independent in the cases studied 
in our paper (see below). Second, we have shown that contributions of gauge modes and 
ghost fields to the full C(0) value do not cancel each other. Hence the reduction of a gauge 
theory to its physical degrees of freedom before quantization only yields the contribution 
of such degrees of freedom to the quantum theory, but is by itself insufficient to describe a 
gauge-invariant quantum theory. Third, we have provided evidence that, when the bound- 
ary 3-geometry does not lead to a well-defined 3-|-l split of the 4-vector potential, some 
inconsistencies occur, i.e. one-loop quantum amplitudes are gauge dependent. 

Interestingly, this seems to complement the Hartle-Hawking programme [5], which 
relies on a Wick-rotated path integral with just one boundary 3-geometry. More precisely, 
on the one hand we know that the use of unitary gauge, 3-1-1 decomposition, Hamiltonian 
formalism and extraction of physical degrees of freedom is necessary to recover the physical 
content of the theory. Moreover, the natural way of implementing the Hartle-Hawking pro- 
gramme [5] involves the consideration of the wave function of the universe in the Lorentzian 
region in terms of physical degrees of freedom (with the subsequent analytic continuation 
to the Euclidean-time region, whenever this is possible). On the other hand, the Euclidean 
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germ from which our universe might originate [5] is the Riemannian 4-manifold which does 
not possess a well-defined 3+1 split [9]. 

Does this mean that it is impossible to carry out the Hartle-Hawking programme in 
a consistent way ? Indeed, it should be emphasized that our present understanding of 
quantum field theory and quantum gravity does not enable one to make a conclusive state- 
ment. Since quantum amplitudes involve differential operators and their eigenvalues, the 
mode-by-mode analysis based on ^-function regularization remains of primary importance. 
Moreover, the analytic continuation back to real, Lorentzian time may be impossible [1]. 
Hence there are cases where we may have to limit ourselves to the elliptic boundary-value 
problems of Riemannian geometry, where one cannot define the notion of time-evolution. 

Other interesting problems remain unsolved as well. In fact, one still has to obtain a 
general proof of gauge invariance in the two-boundary case. Our paper has only focused 
on two particular gauge-averaging functionals. Moreover, one has to repeat our mode-by- 
mode analysis of one-loop amplitudes for spin- 1 fields and gravitation, as well as deal with 
curved background 4-geometries. In the latter case, it is impossible to decouple gauge 
modes without studying fourth-order ordinary differential equations, and one faces the 
technical problem of working out the uniform asymptotic expansions of their solutions. 

The new tools developed in our paper and in the recent literature [1-4,7,9,11] make 
us feel that a complete mode-by-mode analysis of quantized gauge fields and gravitation 
in the presence of boundaries is in sight. Although this would be far from having a good 
theory of quantum gravity, it seems to add evidence in favour of quantum cosmology being 
at the very heart of modern quantum field theory. 
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Appendix 

To clarify our C(0) calculations, some of them are here presented in more detail. We 
first evaluate the contribution to ({0) of non-physical modes in the Lorentz gauge for the 
manifold with one boundary (section 3). For this purpose, let us write the determinant 
(3.4) giving the eigenvalue condition in the explicit form 

^^+1 {{MT+/2) + ^"-2 + j + 1^ C-1 (^^Mt+/2) + + ^"+2) = °- ^^^^ 

The minimal power of M in (Al) is 2n — 1. Thus, to get rid of the fake roots M = we 
should divide our eigenvalue condition by M^""-*^. 

The leading behaviour of /„(nMr) is determined by the exponent e"^"^ multiplied by 
-^y^=- Hence, after taking the logarithm of (Al) divided by M'^^~^ one finds that the only 
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term proportional to logM is (— 2n)logM. Having this expression we can get, by direct 
summation, the result (3.5) for /log. For the calculation of /poie(c>o) it is more convenient 
to re- write (A.l) through derivatives of Bessel functions instead of using the recurrence 
formulae. Thus we have, instead of (Al), the equation 

T+ ( ^^-1 +T'+ \ + T+ ( , w+ \ _n (AO) 

\^ (MT+/2) + )^ (n-l) \^ (MT+/2) + ^"+1 )~^- ^^^^ 

All terms in (A2) have the same leading behaviour determined by exponential functions. 
It is also known that non-exponential terms, after taking logarithms and after expanding 
these in inverse powers of n, are proportional to inverse powers of M and hence vanish in 
the limit M ^ oo [4]. Then equation (A. 2) has the form: exp. terms x ^1 -|- -[^^^jj^ = 0. 
Moreover, it can be shown that also exponential terms do not contribute to /poie(oo) 
[4]. Thus, one should only take the logarithm of ^1 + j and pick out the terms 

2 

which, multiplied by yield terms proportional to 1/n. This leads to the result (3.6) for 

-fpole(oo). 

To calculate /poie(O) it is more convenient to use again the expression (Al) for the 
eigenvalue condition. Taking the limit M — > 0, one finds the limiting form of (Al) as 

1 / I 2\(?i + l)l/l 2 

r(n + 2) Vr(n- 1) ^ r(n)y ^ (n- l)f(^ Vr(n + 1) ^ r(n + 2) 

Thus, inserting the asymptotic Stirling formula for F-functions in the logarithm of the left- 
hand side of equation (A3) and expanding it in inverse powers of n, one finds the result 
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(3.7) for /poie(O). For this purpose, we also use the expansfon of fog(l as a; — > 0, i.e. 

log(l + a;) = Er=i(-ir'^ [!]• 

When one evaluates /poie(oo) and /poie(O) in the Esposito gauge, it is also necessary 

to use the asymptotic expansion oi v s& n ^ oo (cf (5.10) of [2]) 



The contribution of gauge modes to /poie(O) in (2.17) is then obtained by taking the 



coefficient of ^ in the asymptotic expansion as n — > oo of ^ log 
yields the contribution 



(^-l/2)/(^+l/2) ^^ 

r(i/+i/2) 



. This 



9 1 13 111 

I A = \ \ . (^5) 

128 32 360 8 12 48 16 ^ ^ 



Moreover, the contribution of ghost modes to /poie(O) in (2.17) is obtained by taking 

the coefficient of ^ in the asymptotic expansion as n — > oo of ^ log ^ r{u+i) ) then 
multiplying by -2, since the ghost field is fermionic and complex. Hence one finds the 
contribution 

9 11 

Is = — + — - — . (^6) 

128 32 360 ^ ^ 



Thus, the full /poie(O) in (2.17) is given by 



/pole(0)=/A + /B = -^. (^7) 



When one evaluates /poie(oo) and /poie(O) for coupled gauge modes, one finds that only 
K functions at r = t_ and / functions at r = t+ contribute. The corresponding /poie(O) 
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values are obtained as the coefficient of ^ when n ^ oo and M ^ in the asymptotic 

expansions of the terms \ log [^^1/2) ' \ ^'^S t^-^) ' ^ (^+1/2) respectively. 
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